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Synopsis A data reduction procedure for the removal of the effect of geometrical factor from the 
angle dependent specular x-ray reflectivity profile of small samples is described. 

Abstract For small samples, the modification of the XRR profile by the geometrical factors 
manifesting due to profile and size of the beam and the size of the sample is significant. Geometrical 
factors extend till spill over angle which is often greater than critical angle for small samples. To 
separate the geometrical factor, it is necessary to know the spill over angle. Since geometrical factor is 
a smoothly varying function and extends beyond critical angle, it is impossible to determine the spill 
over angle from XRR profile of small samples. We have shown by comparing the normal XRR profile 
of a small sample with the XRR profile taken with a surface contact knife edge on the same sample, 
that the spill over angle can be determined. Thus we have developed a procedure for data reduction 
for small samples and validated it with suitable experiments. Unlike hitherto used methods which 
have drawbacks, this is a self consistent method for data reduction. 

Keywords: Specular x-ray reflectivity; Geometrical factor; Data reduction. 


1. Introduction 

Specular x-ray reflectivity (XRR) is a non destructive technique for characterization of flat 
surfaces or thin films deposited on a flat surface. It mostly determines electron density profile normal 
to the surface which is related to the thickness of different layers and roughness at their respective 
interfaces. Determining such information without physically damaging or altering the sample makes 
the technique attractive for thin film characterization. The theory of XRR is primarily contained in 
Fresnel’s law for reflection of electromagnetic waves at ideal sharp and flat interfaces defined by 
change of refractive index. Layer(s) of thin film(s) on a substrate offer multiple interfaces. The 
reflected x-ray beam from different interfaces, owing to their difference in phase which depends on 
refractive index of different layers and path length, produce Bragg like reflection peaks in the XRR 



profile due to interference effect. Physical roughness, chemical roughness or intermixing between 
layers reduces the contrast at the interfaces due to a relatively gradual variation (transition) of density 
in real interfaces. 

Parratt’s recursive formalism is the procedure to calculate the XRR profile from thin film 
samples by application of Fresnel’s law for sharp interfaces (Parratt, 1954). Novet and Croce model 
incoiporates the effect of roughness by multiplying a factor of exp(-qc>' / 2) to the x-ray reflectivity 
where qz is the vertical wave vector transfer and a is the roughness in angstroms (Nevot and Croce, 
1980). Combination of these two is central and sufficient to describe XRR completely. 

Conventionally, XRR is measured experimentally by a 0-29 scan which records the reflected 
intensity of monochromatic x-ray beam from a thin film sample as a function of 29. The 
experimentally obtained curve is not exactly the same as the Fresnel reflectivity for the following two 
reasons. The first one is due to diffuse scattering which scatters some x-ray photons in the specular 
direction. However the diffuse scattering background is not significant for samples having 
uncorrelated roughness and the method to separate it in the relevant situations is described in literature 
(Romanov et.al., 2010). The second reason is due to the experimental configuration, coined as the 
geometrical factors. It includes modifications in the specular XRR profile due to size of the sample 
and size and profile of the incident x-ray beam. The usual practice is to suitably modify the 
experimentally observed XRR profile by various data reduction procedures described in the literature 
and then to fit it to the theoretical curve calculated using Parratt’s recursive formalism to extract 
numerical values of thickness, roughness and density to different layers. All the available XRR 
analysis softwares generate Fresnel’s reflectivity (modified by roughness factor) to compare to the 
experimentally observed XRR profile. The primary function of the data reduction procedures is to 
remove the effects of geometrical factors and extract the precise Fresnel’s reflectivity from the 
experimentally observed curve. Discussion on whether such a procedure can actually generate the 
precise Fresnel’s reflectivity is lacking in the available literature. The criteria suggested by Gibaud et 
al. to reduce complications in data reduction procedure is to use larger samples and judiciously select 
the size of entrance slit such that: 

Z, > (77 sin ) .(1) 

where L is the length of the sample, T is the aperture of the front slit and 9 C is the critical angle of the 
sample (Gbaud et.al. 1993). When equation (1) does not hold, the XRR profile gets significantly 
modified by geometrical factors which is also known as foot print effect arising due to the size and 
profile of the beam and size of the sample. Reducing the width of entrance slit (T) has limitations. 
(Salah et. al. 2017). At very small apertures, the size of the beam at the sample position becomes 
poorly defined in terms of size and divergence in addition to loss of flux. Furthermore all thin film 
deposition methods cannot produce large samples with desired quality. For example in pulsed laser 
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deposition(PLD), a widely used technique for thin film deposition due to its versatile material 
selectivity and tuneability of deposition conditions, the maximum size with which samples can be 
deposited with desired thickness uniformity is ~ 10mm xlOmm. This makes PLD grown samples 
unsuitable for XRR characterization in the absence of proper data reduction procedure and more so 
when the thin films have smaller value of 9 C . Hence in order to extend the scope of applicability of 
XRR, it is of vital interest to develop a suitable method for correction of geometrical factors for 
situations where it is difficult to satisfy equation (1). It has been proposed in the literature to calculate 
the geometrical factor by calculating the derivative of the reflectivity curve before the critical angle. 
In this work we have argued that this method can be misleading. Traditionally the XRR curve is 
normalized by using the direct beam count or by using an initial guess of the average density of the 
sample. We have discussed in the context of reference (Salah et. al. 2017) and (Gbaud et.al. 1993) 
that both these methods are not correct. 

Even when equation (1) holds, as we have discussed in later sections, the experimental XRR 
profile is not precisely same as Fresnel’s reflectivity. In this paper we have addressed these issues by 
developing procedures to separate geometrical factors from experimental XRR. These procedures in 
general are applicable for samples of all sizes but especially useful for small samples where data 
reduction is critical. Our objective is to develop a data reduction procedure which is applicable to all 
sample sizes and it should not use any of the direct beam count, derivative of the XRR curve, 
physically measured values of sample size or prior knowledge of the optical constants of the thin film 
material to normalize the curve. 

2. Theory 

In this section we elaborate certain important aspects of the experimental geometry of XRR 
which will be used for data reduction. The angular dependence of calculated Fresnel reflectivity will 
exactly correspond to the experimentally obtained reflectivity curve, if the flux available for reflection 
is constant at all angles. The x-ray source provides constant flux with finite and constant beam width 
and a fixed profile. However in the experimental geometry of the 9-29 scan, the angle made by the 
sample of finite size with the beam propagation direction starts from 9=9° with the height of the 
sample so adjusted that it cuts half of the flux from the source. At 9=9° the sample being parallel to 
the beam, does not reflect at all. As 9 increases, the sample surface starts intercepting the beam 
depending on its finite size. This is the reason for a monotonically increasing flux available for 
reflection with increase of angle. The extent of interception of the x-ray beam by the sample is given 
by (Lsin 0) where F is the length of the sample and 9 is the glancing angle of incidence. The 
intercepted beam 0(9) is available for reflection. This can be express as: 

( L/2 )sin6 

®(0)= j P(z)dz .(2) 

(-L/2) sinO 
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where P(z) is the beam profile. If T is the vertical size of the beam then the full beam falls on the 
sample when 9 is such that L sin 9 = T . This angle is called spill over angle 9 S0 . For all 
angles 0 > 9 S0 there is no variation in the flux available for reflection and the experimental reflectivity, 
if normalized correctly, is precisely equal to calculated Fresnel’s reflectivity. For all angles 0 < 9 S0 the 
experimental curve is the Fresnel’s reflectivity modulated by 0(0) given by equation (2). This 
modulation is known as geometrical factor or foot print effect. 


2.1. Geometrical factor due to beam and sample size: 

Fresnel’s reflectivity R(0) of a flat surface is given by: 


W) = 


0-yj0 2 -28 + 2iP 
9 + ^9 2 - 28 + 2if 


(3) 


where 5 and p are the real and imaginary part of refractive index of the sample respectively. 

For simplicity, let us consider that the beam profile P(z) is a rectangular step function of 
width T and height P(z) = 1. With Lsm9 so = T and from equation (1) we get: 


®(<9) = Esin <9, ford<9 so 
= T, for9> 9 so 


(4) 


Equation (4) can be normalized by dividing beam width T to get the geometrical factor g(9). 

g(9) = (L/T)s,m9, for9 < sin -1 (7/ L) 

. .(5) 

= 1, for9 > sin (T / L) 


Where, sin 1 (T/L) = 9 S0 .(6) 

In figure 1 we have plotted the geometrical factor g(9) for different values of (L/T). It can be 
seen in the figure that the spill over angle increases with decreasing (L/T) ratio. Experimental XRR 
profile E(6) can be calculated by multiplying equation (4) to equation (5). In figure 2 we have plotted 
R(9), g(6) and E{9) [=R(#)g(#)J. 2?(<9) is calculated using 8=2.145xlO' 05 and P= 1.262xlO' 06 . It can 

be observed that R(9)=l for 9=0°. In the region of total external reflection (from0 = 0°to 0 = 9°), R(6) 
decreases slowly and monotonically due to absorption factor fi of the sample (Parratt, 1954). This 
feature in not observed in E(6) curve which is the result of modification of R(9) by g(9). The XRR 
measurement produces E(9), and the goal of a data reduction procedure is to produce R(9) from E(9). 
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It can be seen in figure 2 that for 0 > e so , the Fresnel’s reflectivity R(Q) and experimental curve E(6) 
merge with each other. 



Figure 1 Plot of g(0) in log scale calculated for (L/T)=60, 100 and 150. The values of 0 SO in degrees 
are mentioned in the figure. 



0.0 0.5 1.0 

Angle in degrees 


Figure 2 Calculated R(0), g(0) and E(0). 

The actual reflectivity curve contains a part of the direct beam at the values of 0 close to 0°. 
Before starting the 0-20 scan, the usual practice is to align the sample in such a way that it obstructs 
the lower half of the direct beam. Hence during the onset of the scan, the measured intensity as 0 
starts from 0° is the sum of two components. The first one is the unobstructed half beam intensity 
which is being progressively reduced by the increasing obstruction offered by the sample due to 
increase in 0. The second component comes from the reflection of a portion of the initially obstructed 
lower half of the beam by the sample as 0 increases. Furthermore it is well-known that the profile of 
the direct beam obtained from such a scan shows a shift in the angular position of the source due to 
chopping off of the lower half portion of the beam. It is also reported that recording the half of the 
incident flux with the sample parallel to the beam in the half cut positions is also not free form errors. 
Usually the alignment procedures are so critical that the error is often ~10% or larger (Gbaud et.al. 


















1993). Hence we argue that the intensity of the direct beam obtained from the 9-29 scan cannot be 
precisely modelled and parameterized in relation to the actual parameters of the source. For this 
reason we have not used the maximum intensity of direct beam obtained from the experimental curve 
for any of our modelling and subsequent calculations on actual experimental data. Thus we have 
demonstrated that the information contained in the beam which is actually reflected form the sample 
(direct beam is not reflected from the sample) is sufficient to normalize the reflectivity curves. This 
eradicates the effect of alignment errors. 

2.2. Effect of the beam profile on the geometrical factor. 

In this section we have considered that the beam profile of the x-ray source can be a 
rectangular step function, or a Gaussian or a super Gaussian described by equations 7 to 9 respectively 
and shown in the figure 3. The profile functions are: 


T(z)=l, for 
= 0 , for 


\z\<(7 

.( 7 ) 

M>cr 


p ( z ) = Y exp -(z 2 /2cr 2 ) 

P(z)=j exp -(z/V2 a) n 
where n = 4,6,8. 


.( 8 ) 

( 9 ) 


a is the half width at half maximum (HWHM) of the beam and T is the width of the entrance slit. In 
practice, a = T/2 (Gbaud et.al. 1993). 



Figure 3 Rectangular, Gaussian and supper Gaussian beam profiles with order n=4 and n=8. 
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Figure 4 Calculated foot print function g(0) for rectangular, super Gaussian of order 4 and Gaussian 
beam profiles. The sample length and entrance slit width used for calculations are 10,000pm and 
100pm respectively. The values of spill over angle obtained from the calculations are 0.57°, 0.61° and 
0.73° respectively. 

The foot print function g(Q) can be calculated from P(z) as follows: 

(£/2)sin# I (LI 2)sin0 so 

g{&)= | P(z)dz / J P(z)dz .(10) 

-(£/2)sin# / -(L/2)sin0 JO 

It may be noted that equation 10 reduces to equation 5 for rectangular beam profile of width T. Figure 
4 shows the plot of foot print function g(Q) for different beam profiles. The length of the sample and 
the entrance slit width used in the calculation are 10,000 pm and 100 pm respectively. It is observed 
that the slopes of the foot print functions changes according to the type of beam profile. Hence the 
method described by (Baules et.al. 2006) to calculate the sample length from the slope of E(9) below 
the critical angle is only an approximation and more so when we appreciate that this portion of the 
experimental curve is also multiplied with R(6) which is a monotonically decreasing function due to 
finite absorption factor [1 (Baules et.al. 2006). For lower angles, the slopes of the function g(9) remain 
same for all the three type of beam profiles. But at higher angle, they differ. Hence if the critical angle 
0 C of the sample is small, its evaluation by using the foot print function for different beam profiles will 
not vary significantly. For the Gaussian beam profile the curvature of the plot is prominent and 
visible. It can be observed from the figure that the Gaussian beam profile shows highest spill over 
angle (0 SO ). As the beam profile approaches the rectangular box function, 0 SO gradually decreases. 

The awareness that the spill over angle Q so is not only dependent on the physical size of the 
sample and the beam but also on the beam profile, leads to the conjecture that it cannot be calculated 
by simply measuring the physical parameters L and T. As can be seen in figure 4, the shift in spill 
over angle for different beam profiles is as much as 0.16°. Hence it would be incorrect to calculate it 
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from physically measured values of L and T by using equation 6 as suggested in reference (Salah et. 
al. 2017). 

2.3. Classification on the basis of the span of geometrical factor. 



Figure 5 Calculated reflectivity of the same sample with different sample length 

Figure 5 shows calculated reflectivity for two cases. For the curve OA1A2, 0 SO < 0 C and for 
curve OB1B2, 0 SO > 0 C . For reasons described in section 2.1, the curves are normalized to their 
respective maximum instead of the direct beam maximum. We observe that in the curve OA1A2, the 
geometrical factor is present from O to Al and the rest of the curve from Al to A2 is parallel to 
Fresnel’s reflectivity. Since its value at Al is 1 which ideally should have been less than 1 to account 
for finite absorption p in the sample; it can be mathematically represented as A1A2=R(0)/R(0 1 ). 
Available XRR fitting programs calculate R (0) and compare it to the experimentally observed XRR. 
Flowever due care has to be taken when the absorption factor f> is large and 0 SO is very close to but less 
than 0 C . In practice, prior knowledge of 8 and p of the sample are used to generate R(0) and then the 
experimental curve is multiplied with a suitable factor to match with R(0) in such situations. 
Deposition conditions alter density of the samples. There is a fair chance that this change may go 
undetected due to subjective biasing introduced into the fitting procedure by using prior knowledge of 
8 and P to decide the multiplication factor. Flence in such cases, for the sake of theoretical precision, 
we fit the experimental curve to R(0)/R(0i) instead of fitting it to R(0). Here 0i is the angle at which 
experimental reflectivity profile is maximum. 

In the curve OB1B2 the geometrical factor is present beyond Bl. ft is impossible to detect the 
spill over angle. The data reduction procedure developed in this work is categorically useful for this 
situation. However it should be noted that the procedure developed here will be applicable if 0 SO is 
even greater than the maximum range of angle 0s up to which experimental reflectivity signal is 
measurable. 








2.4. 6-20 scan with a surface contact knife edge 


In several commercial reflectometers, beam knife edges are routinely used for reducing the 
beam foot print and reduce air scattering. For small samples it can be used primarily to reduce the 
beam foot print. However due to difficulties involved in the process of alignment it is not easy to 
check the parallelism of the knife edge with the sample and it is also difficult to know the exact 
distance between the knife edge and sample. We use a knife edge made up of normal paper cutting 
blade and make it stand on the sample with its own weight by a simple mechanical arrangement. The 
weight of our knife edge is within 5 grams and is far less heavy than usual commercial knife edges. 

The average distance between the knife edge and the sample is nominal tke ~ . Since the 

nominal sample size is —10,000pm, the (L/T) ratio obtained with this setup is sufficient to yield spill 
over angle far less than critical angle. 



Figure 6 Schematic setup for measurement with a surface-contact beam knife edge 

In the figure 6 we have drawn the schematic diagram of the surface contact knife edge 
arrangement. The allowed beam width can be calculated to be 2 t ke cos 8 . However, taking 9-29 scan 
with this configuration shows that the effect of finite absorption before critical angle is not observed 
in the XRR profile (please refer to Eke( 9) of figure 8). This may be occurring due to small changes in 
the beam size and divergence which can be dependent on the slope error of the sample or small 
misalignments which are unavoidable at these small angles. Due to this distortion occurring before 
critical angle this XRR profile cannot be normalized precisely. But due to high (L/T) ratio, if 
normalized correctly, it will exactly coincide with theoretical reflectivity after critical angle none the 
less. 

2.5. Procedure of data reduction 

In this sub-section we have described a generalized form of data reduction procedure for 
arbitrary beam profile. In section 4 we have elaborated for the specific case when the beam profile is a 
step function. 
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Our data reduction procedure uses two experimental reflectivity profiles of the same sample. 
The first XRR profile is taken without any knife edge and is denoted by E(0). The second XRR profile 
is obtained with a surface contact knife edge denoted by Eke(0). The Fresnel’s reflectivity R(0) is the 
unknown quantity to be evaluated from these information. We assume that E (0) and H K h( 0) can be 


expressed as: 

E(0) = g(0)R(m .(11) 

= 2 .( 12 ) 


Where g{9 ) and g KE (9) are the geometrical factors given by equation (5), Ci and C 2 are the constants 
arising due to normalization. Since spill over angle in Eke( 0) is very small we can consider 
g KF (d) = 1 for all practical purposes. Hence, from equation (12) we get: 

E ke {&) = R(9)C 2 .(13) 

We define E'(6) such that E'(9) = E(9) / (Q / C 2 ); thus from equation (11): 

E'{9) = g(9)R(9)C 2 .(14) 

Writing g(9) explicitly: 

(Z,/2)sin# 

g(&) = Ep f P(z)dz, for9 < sin 9 so 

1 . . z . _ nn 


(-L/2)sin0 .(15) 

= 1, for9 > sin 9 S0 


Where P(z) can be defined by any one of the equations (7),(8) or (9). In order to calculate g(0) it is 
necessary to find the spill over angle. We compare E’(0) and Eke( 0) to find out the spill over angle. 

From equation (15) for ^ - ® so , ~ 1 ; thus by equation (13) and (14); 

E'(9) = E ke (9) for 9 > 9so .(16) 

Equation (16) suggests that ^ ^ and meet each other at ^ = ® so . Thus 0so is determined. 

Satisfying condition that; 

(L/2)sin0 so 

j J p ( z )dz = 0.996.(17) 

(-L/2)sin0 so 

the value of the effective sample length L can be calculated. The above condition guarantees that the 
sample covers 99.6% of the beam at 0so which corresponds to ±3a of the beam profile. (For 
rectangular step function like beam profile, the above condition reduces to (L / T ) sin 9so = 1, from 
where the effective value of (L/T) can be calculated). Now g(0) can be calculated for all values of 0 
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by using the values of L and 0so in equation (15). We identify 9i as the starting value of glancing 
angle and write R(0i) as Ri which is very close to 1. Then we calculate from equation (11) & (15); 

R x C y =E{9 x )lg{9 x ) .(18) 

^ 1 C 2 =R 1 C 1 /(C 1 /C 2 ).(19) 

Using equation (11) & (18) we calculate; 

T(0) = R(d)/R x = E(0) I [g(6)R x C x \ .(20) 

And using equation (12) & (19) we calculate; 

T KE (6) = R(0)/Rl = E ke (6) / [R1C 2 ].(21) 

Since Ri=R(0i) ~ 1, both ¥(0) and T , K e(0) can be used independently used for fitting with Fresnel’s 
reflectivity. For the sake of generality and theoretical accuracy, we fit T (0) and 1 P K k (9) with R(0)/Ri. 

It should be noted that if L and T are the physically measured dimensions of the sample and 
the beam respectively, then the maximum value of 0so considering a Gaussian beam profile can be 
calculated by the following formula for a quick approximation. 

£“ ax =shT 1 (1.2677 Z).(23) 

3. Experimental 

3.1. Sample Preparation 

We have deposited E 1 - 2 O 3 thin film samples on GaAs substrate by pulsed laser deposition 
(PLD) technique using Excimer laser of wavelength 248nm. The dimension of the substrates was kept 
at lOmmXIOmm. It’s a well know fact that desirable thickness uniformity is not achieved by PLD if 
bigger substrates are used for deposition. Five samples were deposited by varying the substrate 
temperature and ambient oxygen pressure during deposition. Laser energy, repetition rate and 
deposition time were similar for all the samples. The details of the samples are listed in table 1. 

Table 1 Deposition environment of various samples. 


Sample Name 

Substrate temperature 
in oC 

Chamber pressure in 
mbar 

S 1 

550 

9x10-5 

S 2 

600 

9x10-5 

S 3 

650 

9x10-5 

S 4 

550 

5x10-3 

S 5 

550 

5x10-2 
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3.2. Measurements 


Normal XRD measurements in 9-29 geometry were done for all the samples. It was found that 
the samples were crystalline and variation of deposition temperature has no visible effect on the 
crystal structure of the material till ambient oxygen pressure was maintained below 5x19 3 mbar. 
However when oxygen partial pressure was increased to 5x19 2 mbar, the (222) peak of Er 2 03 shifted 
to lower 29 value which corresponds to the actual lattice parameter of bulk Er203. This suggests that 
the lattice contracted for lower oxygen partial pressure. 



Figure 7 XRD pattern of the sample deposited at 5xl9' 3 mbar pressure is plotted in black and that of 
the sample deposited at 5xl9' 2 mbar is plotted in red. 

Figuer 7 shows the change of lattice parameter of Er203 layer in sample 4 and 5, due to 
change in ambient oxygen pressure. This change in lattice parameter must reflect in change in density 
of the samples which should be detectable in the XRR measurements. With the method of data 
reduction described in section 2, we carried out the XRR analysis of all the samples. The hard x-ray 
characterizations of all the samples were done using the grazing incidence x-ray reflectivity (XRR) 
technique on a Bruker Discover D8 reflectometer and diffractometer with Cu Ka radiation (A. = 1.54 
A) and using a Kp filter. The 9 and 29 values were aligned to better than 9.992°. Specular reflectivity 
was measured in 9-29 geometry over the range of 9° to 4° with a step size of 9.995° for all the 
samples. As described in section 2.4, we also recorded the XRR measurements on the samples using 
surface contact knife edge. 

A program was developed for data reduction in MATLAB. The processed data were then 
used for reflectivity fitting using Parratt32 reflectivity fitting software. 

4. Analysis 

The direct beam part of the normal XRR profile and the surface contact knife edge XRR 
profiles were removed and then they were normalized with respect to their respective maximum. They 
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are plotted in figure 8 and demarcated by E{0) and Eke (i 0 ) respectively. We assume a step function 
like beam profile described in equation 7. 



Figure 8 Reflectivity measured with and without Kn ife edge plotted against angle. 



Figure 9 E(0) (black curve) is divided by a constant to be identified as (C1/C2) so that the high 
angle part of both the curves coinside 

Referring to figure 8, we identify the x-coordinate of the first data point of F(0) as 0 X . The y- 
coordinate of the first data point of F(0 ), by equation 11 is: 

E(0 \) = g(9\ )RA (23) 

Here we identify the Freshnel reflectivity at 0 X to be K } . 

As discussed in section 2, E ke (0) is same as R(<9)for higher angle. We divide E (<9) by a 
suitable factor to obtain E \0) so that Eke (0) and E\0) coincide at higher angle. We identify the factor 
to be (C1/C2). In figure 9 we have plotted both the curves after performing this operation. The angle 
at which both these curves meet is identified as the spillover angle. 
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We apply (L / T) sin 9 so = 1; to calculate (L/T). g(0) is calculated from equation (5) by using the 

value of (L/T). Thus the value of g(0i) is found which is used in equation (18) to find out the value of 
R1C1. Then R1C2 can be found by equation (19). 

T(6() and T^)#) are determined by using equations 18 to 21and the values of g(0), R1C1 
and R2C2. We have plotted E(0) and 'T(O) in figure 10. It can be observed that ¥(0) which is derived 
from E(9) shows the monotonically decreasing trend which is the effect of finite absorption before the 
critical angle. 

It is worthwhile to study the behaviour of the function g(9)R l C l with respect to the spill over 
angle 9 so . Change in 9 so changes g(9) through equation 5 and 6. So, following equation 5, we write: 

g(9)R l C l =E(9 l )sin9/sin9 l , for 9<9 so 

.(25) 

= m> M 0>9 SO 

Let us consider the case: 9 < 9 S0 . In the RHS of the above equation, E(9 { ) and 9 X are constants taken 
from the experimental curve. They are independent of 9 S0 . So the REIS is independent of spiller over 
angle. But in the LHS, g{9) is calculated from 6 S0 . If there is any uncertainty in evaluation of the spill 
over angle, there will be uncertainty in the value ofg(<9). In such situation, to satisfy LHS=RHS, the 
factor RiCi in the LHS will take an appropriate value. In other words g(9) and RiCi conspire to keep 
the product g(9)R x C x a constant independent of spill over angle. It is this feature of the procedure that 
makes it appropriate for data reduction even when the spill over angle is greater than the maximum 
angle up to which data can be measured (i.e. for very small samples sizes). The function g(9)R l C 1 , for 
different values of spill over angle are plotted in figure 10. It can be observed that this function is 
unique for 0 < 0 SO for all values of 0 S o. 
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Figure 10 E(9) and l P(0) are plotted. The curves of g(9)RlCl for different values of spill over angle 
are also plotted. The spill over angle used for calculation of ¥(9) is 0.88 degrees. 



Figure 11 The normalized reflectivity curves ¥(9) and V FKE(9) are plotted. 

Figure 11 shows the final normalized curves after data reduction. The constant C 2 is related to 
the spill over angle through equations 6, 18 and 19. Both these curves can be fitted with the Fresnel 
reflectivity to produce identical results. The mismatch of these curves before critical angle shows that 
the knife edge modifies the XRR profile before critical angle. 

5. Results and discussion 

Reduced data for all the five samples were taken for fitting in Parratt32 reflectivity fitting 
software. The XRR with the knife edge show clear thickness oscillations where as these are smeared 
out in the XRR profile without knife edge. The reason for this is the increased resolution of the 
measurement due to reduction of beam width by the knife edge. We have confirmed by fitting both 
these data that both of them correspond to same density variation model with different experimental 
resolution. The fitting are shown in figure 12. In table 2, we have mentioned the fitting parameters 
applicable for both curves. 


Table 2 XRR fitting parameters 



Thickness 

In angstroms 

Real part 

of p 

Imaginary 

part of p 

Roughness 

in angstroms 

Interface 

93.13 

5.5E-6 

1.6E-7 

7.6 

Er 2 0 3 

643.6 

5.8E-5 

4.7E-6 

6.7 

Interface 

61.04 

3.9E-5 

5.2E-6 

25.5 

substrate 

Inf 

3.8E-5 

1.1E-6 

4.1 
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Figure 12 Fitting of *F(0) and 'FKE(0). Both of them fit with same density variation model but 
with different experimental resolution. So, both of them represent same micro structural information. 

Sample 5 was deposited at 5x1 O' 2 mbar pressure and the XRD measurements shown in figure 
7 suggest that there was an increase of lattice parameter for this sample. The O-scan in XRD 
measurements ( graph not shown here) confirm that the film was not epitaxial. Hence it has no 
relation to the lattice parameters of the substrate which indicates there is no biaxial strain in the 
sample. So this change in lattice parameter must correspond to a change in density of the sample. The 
scattering length density from XRR measurements shown in figure 13 also show that the density of 
this sample has decreased by ~6% than the other four samples. Furthermore it is worthwhile to 
mansion that the density of samples SI to S4 are ~4.6% higher than the actual bulk density of FiyOi 
where as density of sample 5 is 1.6% lower. The detection of this change of density would not have 
been possible by the prevailing method of data normalization where the experimenter uses the 
available value of 8 and p to normalize the XRR profiles and thus introduces an unwanted subjective 
bias into the fitting procedure. The small change of density due to variation of deposition parameters 
would have gone unnoticed. 
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Figure 13 Scattering length Density variation modes obtained from fitting of XRR data of the 
Er203samples. 

We have also performed the data reduction assuming the beam profiles to be a Gaussian or 
super Gaussian described by equations 8 and 9. They yield similar result. The difference in the 
deduced values of 8 and p were negligible. This is as a result of insignificant difference between g(0) 
values calculated for these beam profiles at lower angle as is evident in figure 4. 

6. Conclusion 

We have proposed a method of data reduction for the XRR measurement of small samples. 
We have not used direct beam count, derivative of the XRR curve, physically measured values of 
sample size or prior knowledge of the optical constants of the thin film material to normalize the 
curve. We have elucidated the effects of size and profile of the x-ray beam and size of the sample on 
the geometrical factor. Using this method, geometrical factors could be successfully removed from the 
data to yield Fresnel reflectivity. Determining spill over angle was the key issue in determining the 
geometrical factor. The modus operandi described in this work for determining spill over angle by 
comparing usual experimental XRR profiles with the XRR profile taken with a surface contact knife 
edge is a heuristic solution that is well explained by intuitive modelling of the geometrical factor. This 
is the main aspect which makes it different and effective from the available literature. We have shown 
that theoretical Freshnel reflectivity R(9) can never be extracted from the experimental reflectivity 
curve. Only R(0) /R(0i) can be extracted by using a legitimate data reduction procedure. Flence we 
propose to do away with the prevalent procedure of fitting the experimental reflectivity with R(9) in 
the favour of fitting it with R(9)/R(9i), where 9i is the lowest angle after removing the direct beam 
part. The success of the method is demonstrated on Er 2 0 3 thin film samples. Though it will be 
incorrect to draw any precise numerical connection (it is beyond the scope of this work) between the 
variation of lattice parameter found from XRD and the density of the material found from XRR, it is 
none the less the success of the method to detect such precise density variation when there is a change 
of lattice parameter. 
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